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This  final  report  for  AFURL  eontraet  FltM>28-7.‘L(’-02Hfi  constitutes  a small  increment 
m a continuing  series  of  surface  wave  studies  of  the  Farth's  crust  and  upper  mantle  at  the 
Universitv  of  Texas  at  1 (alias.  These  studies  have  been  based  upon  the  use  of  modern  methods 
for  t lie  analysis  of  dispersed  clastic  waves;  a significant  proportion  of  these  developments  have 
originated  within  our  research  jjroup.  The  present  report  deals  primarily  with  surface  wave 
studies  at  NORSAR.  which  were  performed  under  the  auspices  of  the  Air  Force  Cambridge 
Research  Laboratories.  Various  portions  of  this  material  were  reported  during  and  as  part  of 
an  mv  ited  lecture  “Analysis  of  ( )seil  hit  ions”  ( Landisman  l!*7f>(.  to  the  X VI  ( leneral  Assembly 
of  the  International  I’nion  of  (ieodesy  and  Geophysics,  which  was  convened  in  Grenoble, 
Fram  e,  during  Au^rnst-September.  ll(7fi.  Additional  material  in  this  rep*  rt.  and  o* her  related 
topics  will  appear  in  other  studies,  such  as  that  by  Nyman  and  Landisman  ( 1 i*7<»). 
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Final  Technical  Report 


Introduction 


Recent  research,  including  seisnidlojricjil  ami  gravity  studies  of  t.ieOsloOraben  (Rambert 
and  Smithson.  l!>7ol  and.  more  particularly  of  the  NORSAR  site  (Berteussen.  l‘F5)  have 
tended  to  emphasize!  he  tfeoloybe  and  tectonic  complexity  of  the  region  near  thisarray  site,  Aki. 
et  al  l llBn).  for  example,  have  presented  evidence  for  pipe-like  tfeo  logical  structures  under 
NORSAR.  which  were  interpreted  to  he  the  result  of  plastic  deformation  of  the  local 
hthosphere.  The  discussion  of  these  anomalies  included  consideration  of  possible  consequences 
of  the  a>  *ent  of  marina  associated  with  Permian  volcanism  in  this  portion  of  the  Oslo  (Jrahen. 
The  widely  observed  inhoniotfeneities  in  wave  propagation  for  the  compressional  waves 
recorded  at  NORSAR  are  related  to  inhomo^eneities concentrated  in  the  lithosphere  beneath 
the  array  ( Berteussen.  I!»7f,l.  The  scale  sizes  for  the  inferred  lithospheric  inclusions,  as  judged 
from  spectral  analyses  of  these  perturbations  of  the  compressional  wave  field,  ranue  from 
several  tens  to  over  a hundred  kilometers  (ibid). 

The  propagation  of  lon^r  period . tel esei stoic  Rayleigh  and  I ,ove  waves  tends  to  be  relatively 
indi-turhcd  by  lithospheric  variations  near  the  smaller  end  of  this  size  distribution,  since  the 
smallest  lengths  of  these  surface  waves  certainly  exceed  some  several  tens  of  kilometers.  The 
otrm-r  elements  of  this  distribution  of  in  homogeneities,  however,  appear  to  represent  an 
extremely  serious  problem  for  propagating  elastic  surface  waves.  As  a result,  lateral 
refraction  and  reflection,  scattering  offprreat-cirde  and  non-least-time  arrivals  seem  to  he  at 

least  as  common  for  teleseismic  surface  wave  arrivals  at  NORSAR  as  they  are  for  teleseismic 
body  w ave  arrivals. 


Preceding  page  blank 


-5- 


1 Re\  1CU 


t 

■ 


i 


j 


Our  NOKSAK  surface  wave  studies  have  centered  on  several  main  types  of  problems. 
F irst  w e concentrated  on  the  determination  of  and  eompensat  ion  for  long  period  instrumental 
effects  I hose  efforts  received  direct  assistance  from  laboratory  determinations  of  the 
magnification  curve  ((indzm.  ldbHl  and  synoptic  measurements  of  the  appropriate 
instrumental  constants  at  the  NORSAR  lonjr  period  site,  kindly  supplied  by  ,1.  YVoolson  of  the 


1 idedyne-t reotech  Seismic  Data  Laboratories  m Alexandria.  Compensation  for  the 
magnification  and  phase  shift  according  to  the  frequency  domain  equalization  techniques 
proposed  by  1 lerejrewsk)  ( 1!)71).  was  jrreatly  favored  by  the  fact  that  the  data  were  originally 
recorded  as  digital  time*  series.  1 hey  therelor  possess  a far  greater  dynamic  ranjre  than  is 
ordinarily  available  for  older,  visible  recordings. 

Frequency  domain  compensation  for  the  recording  system  response  (ibid)  broadens  the 
usable  bandwidth  of  recorded  signals  io  a remarkable  decree.  It  also  distinctly  extends  the 
ranye  of  applicability  of  these  surface  wave  studies.  For  example,  whereas  the  original 


recording  appear  to  consist  almost  entirely  of  signals  having  periods  of  from  twenty  to  thirty 
seconds  because  the  system  magnification  is  sharply  peaked  near  periods  of  twenty  five 
seconds,  the  compensated  recordings  have  permitted  the  study  of  dispersion  data  from  periods 
ranvrinvr  from  less  than  fifteen  seconds  to  more  than  two  hundred  seconds. 

Further  efforts  were  concerned  with  determination  of  thedispersionof  Ray  leitfh  waves  for 
the  hurasian  land  mass  for  paths  close  to  potential  test  sites.  One  of  tin*  most  suecessfu  1 ways  to 
accomplish  this  task  involves  the  use  of  frequency  'one  analysis  (FTAN)  of  the  time  series 
1 1 .and istnan.  ot  al..  l!Mi!b  recorded  at  tin*  various  NORSAK  lonyr  period  sites,  after 
instrumental  compensation,  as  discussed  above.  FTAN  processing  permits  determination  of 
the  observed  average  ffroup  velocity  for  the  path  between  the  source  and  the  array  site. 


. 
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A mtii's  of  publications  has  described  the  multiple  filter  process,  which  often  has  been 
used  for  frequency-time  analysis (Dziewonski.  Bloch,  and  Landisman.  1969;  Landisman,  etal.. 
19*>9:  Dziewonski.  Landisman.  et  al..  1968).  This  process  yields  displays  of  “instantaneous” 
spectral  amplitude  and  phase  as  functions  of  log  period  and  linear  group  velocity.  During  this 
phase  of  our  research.  Dr.  Douglas  ('.  Nyman  became  impressed  with  the  fact  that  equalization 
of  the  velocity  and  period  resolution  of  such  displays  would  be  a distinct  improvement  over 
previous  analyses  based  upon  a constant  relative  band  width,  a . (Dziewonski.  Bloch,  and 
Landisman,  19*19).  The  display jequaMzcd  multiple  filter  proeess  (Nyman  and  Landisman, 
197*1)  employs  a locally  variable  bandwidth 


*(Tj,vk)=/.  fv  Tj/  ( 4r  t vkI 2), 


where 


and 


L = epicentra'  distance, 

• v = linear  step  in  analysed  group  velocity  between  adjacent  rows  of  display, 
period  being  analysed, 

fractional  change  in  period  between  adjacent  columns  of  display, 
group  velocity  being  analysed. 


tj 


I he  display  equalized  filler  thus  defined  is  a function  of  the  particular  period  or  velocity  of 

interest  in  contrast  to  the  previously  used  constant  relative  bandwidth  filter  function,  a.  which 
tended  to  have  excessively  great  velocity  resolution  at  short  periods  and  low  group  velocities, 
and  excessive  frequency  resolution  at  long  periods  and  high  group  velocity  values.  Compare 
1-  igures  la  and  lb.  The  resolution  is  also  independent  of  the  shapeof  thedispersion  curve  and  of 
the  spectral  excitation  (Nyman  and  Landisman,  1976).  The  use  of  this  form  of  frequency-time- 
analysis  has  enabled  us  to  make  accurate  estimates  of  the  observed  Rayleigh  wave  disprison  at 
each  long  period  site,  and  has  permitted  the  construction  of  time  variable  filters  (Landisman. 
et  al..  1969;  Dziewonski.  Landisman.  et  al„  1968)  for  isolation  of  the  fundamental  Rayleigh 
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mod*-  signal  on  the  vertical  and  longitudinal  components  at  each  NORSAK  long  period  site 
(Figure  h-i.  Further  purification  of  these  signals  wasderived  from  the  methods  based  upon  the 
use  of  three  component  recordings  to  emphasize  those  signals  arriving  from  the  proper 
azimuth,  with  the  expected  retrograde  elliptical  motion,  and  with  expected  amplitude  ratios 
between  the  longitudinal  and  vertical  components. 

All  of  these  criteria  plus  the  coherency  of  the  longitudinal  and  vertical  arrivals  were  used 
to  improve  the  purity  of  the  dispersed  Rayleigh  mode  arrhals.  especially  the  purity  of  their 
Phase  functions.  We  previously  had  attempted  todetermino  the  phase  velocity  d ispersion  curve 
in  the  vicinity  of  the  NORSAK  array  site,  through  the  analysis  of  “raw”  observations  of  the 
dispersed  wavetrams.  I he  results  for  both  ordinary  and  high  resolution  frequency 
wavenumber  processing  (MeCowan  and  Liritz.  1W8)  of  such  data  were  all  but  meaningless. 
Careful  "tud\  of  the  long  period  horizontal  recordings  showed  that  pseudo-random  lateral 
refractions  are  an  important  component  of  the  problem:  two  dimensional  frequency  wave- 
number  analyses  showed  that  azimuthal  errors  of  the*  order  of  twenty  degrees  can  be  observed 
e\on  for  time  variable  filtered  data,  and  for  periods  that  exceed  100  seconds.  This  is 
surely  a strong  confirmation  of  body  wave  studies  (Aki.  et  a]..  1075)  as  well  as  of  regional 
geological  and  geophysical  investigations  (Rami, erg  and  Smithson.  1075).  that  lateral 
in  homogeneities  having  scale  sizes  jn  excess  of  one  hundred  kilometers  can  significantly 
influence  wave  propagation  in  the  vicinity  of  the  NORSAR  array 

These  and  similar  effects  naturally  become  still  more  important  at  shorter  periods  and 
shorter  wavelengths.  Monitoring  of  the  longitudinal-vertical  coherency  and  of  the  azimuthal 
deviations  of  the  dispersed  arrivals  appears  to  provide  a reasonable  means  of  predicting  those 
portions  of  the  spectrum  mat  correspond  primarily  to  the  great-eirele-azimuth.  least-time 
dispersed  arrivals.  These  predictions  tend  to  be  confirmed  in  most  cases  by  the  two 
dimensional  frequency  wave  number  analyses  of  the  horizontal  recordings.  Viz:  those  portions 
of  a dispersed  wav, drain  that  arrive  from  azimuths  close  to  the  great  circle  direction  of 
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approach  also  exhibit  good  longitudinal-vertical  coherency  and  expected  values  of  ellipticrty. 
The  converse  statement  is  usually  observed  as  well.  This  is  of  distinct  significance,  since  the 
relatively  rapid  analysis  of  a few  sample  traces  can  thus  be  used  with  considerable  confidence 
as  a guide  to  the  use  of  slower,  more  expensive  frequency  wavenumber  analysis  of  data  from 
the  entire  array. 

The  improved,  though  still  fairly  rough  phase  velocities  defined  by  this  process  should  be 
able  to  be  refined  by  iterative  analyses  of  a large  suite  of  events  arriving  at  the  array  site  from  a 
w ide  range  of  azimuths.  This  investigation  will  be  continued  in  future  studies,  depending  in 
part  upon  the  availability  of  funds. 
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PERIOD 


FIGURE  1A 

F r e q u e n c y - 1 i m e analysis,  u s i n k constant  relative  bandwidth 
for  NORSAR  recording  of  Java  event,  mh=6.1,  1 May  1%1,  a=98°. 


uuntm 


I*  roquency-time  analysis  of  sami*  data  as  in  Figure  la,  processed  by  display 
equalized  filter. 


PERIOD 


FIGURK  1C 

Frequency-time  analysis  as  in  Figure  lb.  after  time  variable  filtering  to 
extract  dispersion  for  the  fundamental  Rayleigh  mode. 


